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ABSTRACT 

We present a study of the variation of spatial structure of stellar populations 
within dwarf galaxies as a function of the population age. We use deep Hubble Space 
Telescope/ Advanced Camera for Surveys imaging of nearby dwarf galaxies in order to 
resolve individual stars and create composite colour-magnitude diagrams (CMDs) for 
each galaxy. Using the obtained CMDs, we select Blue Helium Burning stars (BHeBs), 
which can be unambiguously age-dated by comparing the absolute magnitude of indi- 
vidual stars with stellar isochrones. Additionally, we select a very young (< 10 Myr) 
population of OB stars for a subset of the galaxies based on the tip of the young 
main-sequence. By selecting stars in different age ranges we can then study how the 
spatial distribution of these stars evolves with time. We find, in agreement with pre- 
vious studies, that stars are born within galaxies with a high degree of substructure 
which is made up of a continuous distribution of clusters, groups and associations from 
parsec to hundreds of parsec scales. These structures disperse on timescales of tens to 
hundreds of Myr, which we quantify using the two-point correlation function and the 
Q-parameter developed by Cartwright & Whitworth (2004). On galactic scales, we can 
place lower limits on the time it takes to remove the original structure (i.e., structure 
survives for at least this long), t cvo , which varies between ~ 100 Myr (NGC 2366) 
and ~ 350 Myr (DDO 165). This is similar to what we have found previously for the 
SMC (~ 80 Myr) and the LMC (~ 175 Myr). We do not find any strong correlations 
between t cvo and the luminosity of the host galaxy. 

Key words: galaxies - dwarf; galaxies - individual (DDO 165, Holmbcrg II, IC 2574, 
NGC 784, NGC 2366, NGC 4068) 



1 INTRODUCTION 

Young stars are observed to be born in a clustered, hierar- 
chical or fractal distribution within galaxies, and the goal of 
the present work is to measure the timescale over which 
stellar structures in dwarf galaxies disperse. In Gieles et 
al. (2008; hereafter G08) and Bastian et al. (2009; here- 
after B09) we investigated the stellar distribution in the 
SMC and LMC, respectively. In these works we developed 
a series of statistical techniques with which structures can 
be quantified and their evolution traced as a function of 



time. We found that for the SMC and LMC stellar substruc- 
ture was erased gradually, being indistinguishable from the 
background distribution by 75 Myr and 175 Myr, respec- 
tively. These timescales are comparable to the crossing time 
(tcross = i?gaiaxy/o" where a is the stellar velocity disperison) 
of the galaxy which led to the conclusion that general galac- 
tic dynamics were responsible for the removal of structure, 
rather than local effects such as popping star clusters (also 
known as infant mortality, i.e., cluster expansion due to the 
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loss of residual gas left over from the star-formation process 
- Boily & Kroupa 2003; Bastian & Goodwin 2006). 

One of the major limitations of the previous works was 
in defining samples of nearly coeval stars (i.e. stars of simi- 
lar age). Previously, we selected groups of stars from regions 
in the colour-magnitude diagram along the main-sequence 
(MS) and then carried out simulations to estimate the mean 
age and spread of stars within that region of space in the 
colour-magnitude diagram (CMD). This resulted in a nearly 
coeval population only for the youngest ages, whereas re- 
gions further down the CMD were made up of older stars 
as well as some contribution of younger stars. Our results 
were tested against the evolution of the distribution of stel- 
lar clusters with known ages, and overall we found good 
agreement. 

The post-MS evolutionary phase of young (<1 Gyr) 
intermediate mass (> 3 Mq) stars provides an alternative 
method of age dating that avoids the age degeneracies of MS 
stars. Of particular importance are core blue helium burning 
stars (BHeBs) which have a one-to-one relationship between 
age and luminosity (e.g., Bertelli et al. 1994), making it pos- 
sible to assign ages to individual BHeBs on an optical CMD. 
The short lifetimes of BHeBs, ~ 10% of their MS lifetimes 
(e.g., Chiosi, Bertelli, & Bressan 1992), allow for a clear iden- 
tification of coeval populations with minimal confusion be- 
tween stars of different ages (Dohm-Palmer et al. 1997). For 
a given age, BHeBs are ~ 2 magnitudes brighter compared 
to MS stars, thus providing higher photometric precision and 
longer look back times (~ 500 Myr; e.g., Dohm-Palmer et 
al. 1997, Weisz et al. 2008). Although there are fewer BHeBs 
than MS stars, the precise age-dating, possible with BHeBs, 
compensates for the reduction in numbers and provides a 
more secure final sample of stars. In those cases where the 
photometry is sufficiently deep, the star formation histories 
derived from the BHeB stars are consistent with those de- 
rived from the main sequence stars (e.g., Dohm-Palmer et 
al. 2002). 

Imaging of resolved stellar populations, specifically 
BHeBs, in nearby star forming dwarf galaxies, i.e., dwarf 
irregular galaxies (dlrrs), could reveal fundamental pro- 
cesses driving the evolution of stellar structures. Unlike 
more massive galaxies, dlrrs are close to solid body rotators 
(e.g., Skillman et al. 1988, Skillman 1996), implying that the 
locations of stars are not significantly disturbed by external 
mechanisms such as sheer or spiral density waves. Combined 
with precise ages of BHeBs, the dynamic simplicity of dlrrs 
may provide new insight into the intrinsic evolution of stellar 
structures. 

Ideally, one would like to study not only the timescale 
over which stellar structures are erased, but also how this de- 
pends on spatial scale. Sub-cluster (< 1 pc) structures are 
erased quickly (Goodwin & Cartwright 2004) while larger 
structures take more time to disperse. The exact dependence 
of the structural evolutionary timescale and the spatial scale 
probed will be essential in confirming the cause of the evo- 
lution. 

There have been many recent studies on quantifying 
the distribution of point-like (i.e., stellar) structures on a 
variety of physical scales. Near- and mid-IR studies of star- 
forming regions on scales of 0.1 — 10 pc have revealed a wealth 
of structural information (Lada & Lada 2003; Cartwright 
& Whitworth 2004; Gutermuth et al. 2005, 2008, 2009; 



Schmeja & Klessen 2006; Koenig et al. 2008, Kraus & Hille- 
brand 2008) with many regions appearing centrally concen- 
trated and relaxed while others are still largely fragmentary. 
On intermediate spatial scales (1 — lOOpc) extragalactic stud- 
ies on objects like NGC 346 in the SMC have revealed that 
star-forming regions are made up of many individual clus- 
ters and also distributed star-formation (e.g., Sabbi et al. 
2007; Schmeja, Gouliermis, & Klessen 2009). On galactic 
scales, it appears that there is not a characteristic size for 
star- forming regions, and that there exists a continuous hi- 
erarchy of structure from pc to kpc scales (Elmegreen & 
Elmegreen 2001, Elmegreen et al. 2006; Bastian et al. 2005, 
2007, 2009; G08; Odekon 2008). The framework that is de- 
veloping is that stars form in a highly structured distribu- 
tion on all scales probed so far, and that these structures 
are erased on approximately the crossing time, regardless of 
scale. 

The number of nearby dwarf galaxies with high quality 
optical imaging available has significantly increased in the 
past decade, thanks to a number of Hubble Space Telescope 
survey programs (e.g., Tully et al. 2006, Weisz et al. 2008, 
Dalcanton et al. 2009). The excellent spatial resolution af- 
forded by HST allows for precise photometry and age dat- 
ing of individual BHeBs. In this paper, we have selected six 
star forming dwarf galaxies in the Local Volume (Weisz et 
al. 2008, McQuinn et al. 2010) in which we will explore the 
evolution of stellar structure using a statistically significant 
number of BHeBs. 

This paper is organised as follows. In §[2] we introduce 
the observations used in the present study. The technique 
of age-dating BHeB stars is discussed in § [3] and in § [4] we 
introduce the methods used to quantify the evolution of spa- 
tial structures along with their caveats. The results for each 
individual galaxy are given in § [5] and in § [6] we summarise 
and discuss our main findings. 



2 OBSERVATIONS AND PHOTOMETRY 

Observations of all galaxies were taken with the Advanced 
Camera for Surveys (ACS, Ford et al. 1998) aboard HST. 
DDO 165, Holmberg II, NGC 2366, and IC 2574 were 
observed in the F555W and F814W filters as part of a 
survey of the M81 Group dwarf galaxies (GO-10605 - 
Weisz et al. 2008; also GO-9755, PI Walter, was used 
for IC 2574), while NGC 4068 (GO-9771 - Karachentsev 
& Kashibadze 2006) and NGC 784 (GO-10210 - Tully et 
al. 2006) were both observed in the F606W and F814W fil- 
ters. The data we use in this paper have been processed as 
described in Weisz et al. (2008) and McQuinn et al. (2010). 
Here, we briefly summarize the data reduction process. 

Images were all downloaded from the HST archive fol- 
lowing standard reduction and processing with the HST 
pipeline. To ensure uniformity among our data, photome- 
try of all images was performed using dolphotEI, a ver- 
sion of HSTPHOT optimized for the ACS (Dolphin 2000). 
Cosmic rays, hot pixel residuals, and extended objects were 
all rejected based on their brightness profiles. The remain- 
ing objects were subjected to further photometric cuts. 

1 http:/ /purcell. as. arizona.edu/dolphot/ 
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Specifically, we considered well measured stars to have a 
signal-to-noise ratio > 5, \sharpi + sharps] < 0.27, and 
crowdi + crowds < 1.0. The definitions of these photomet- 
ric criteria can be found in Dolphin (2000) and Dalcanton et 
al. (2009). Full descriptions of data reduction for NGC 4068 
and NGC 784 can be found in McQuinn et al. (2010), while 
identical details for DDO 165, Holmberg II, NGC 2366, and 
IC 2574 are presented in Weisz et al. (2008). 

Galaxies considered in this paper are a subset of ~ 25 
star forming dwarf galaxies presented in Weisz et al. (2008) 
and McQuinn et al. (2010). From this larger sample, we first 
eliminated galaxies whose CMDs did not have a sufficient 
population of BHeBs. From the remaining sub-sample, we 
then considered only galaxies which have clearly separable 
MS and BHeB sequences. Simulated CMDs (Gallart, Zoc- 
cali, & Aparicio 2005, Weisz et al. 2008) show that the MS 
and BHeBs are typically separated by several tenths of a 
magnitude in colour, enough to trivially identify them as 
two separate sequences. However, in the event of intense pe- 
riods of SF or 'starbursts', a feature which describes many 
of the galaxies in the initial sample of ~ 25 (McQuinn et 
al. 2010), internal extinction effects can smear out the sepa- 
ration between the MS and BHeBs (e.g., Cannon et al. 2003, 
McQuinn et al. 2009, 2010). Additionally, higher values of 
metallicity lead to higher dust-to-gas ratios with greater con- 
comitant differential extinction. This effect can make it dif- 
ficult to discern the difference between individual young MS 
and BHeB stars. As discussed in § 14.31 and Appendix B in 
B09, inclusion of young MS stars in the BHeB sample can 
introduce false structure detections, potentially biasing the 
results. Thus, we eliminated galaxies without clearly sepa- 
rated MS and BHeB sequences. Finally, we attempted to 
include larger, actively star-forming galaxies to our sample, 
in particular NGC 4449 and NGC 5253. Unfortunately, these 
had to be eliminated due to completeness issues in the high 
surface brightness central regions of the galaxies for the older 
(fainter) BHeB stars. Application of these criteria to the ini- 
tial sample resulted in the six galaxies we consider for this 
study. 

The CMD of Holmerg II is shown in Fig. [T] and its main 
features are indicated. Additionally, the observed CMDs of 
all galaxies used in the present student are shown in Fig. [2] 



3 AGE DATING METHOD 

Precise age dating is essential to understanding the evolu- 
tion of stellar structures. Fortunately, luminous BHeBs have 
high degrees of photometric accuracy. Combined with the 
one-to-one age luminosity relationship predicted by stellar 
models (e.g., Bertelli et al. 1994), it is possible to assign 
unambiguous ages to BHeBs < 500 Myr old (e.g., Dohm- 
Palmer et al. 1997). In this paper, we use the stellar models 
of Marigo et al. (2008), which include the massive star mod- 
els of Bertelli et al. (1994) transformed into ACS instrumen- 
tal magnitudes. 

For our choice of observational wavelengths, the BHeB 
sequence is a primarily vertical sequence of stars located red- 
ward of the young MS, which generally occupies a relatively 
unique position on the CMD (see Fig.QQ. The exact location 
and shape of the BHeB sequence depends upon both metal- 
licity and differential foreground and/or internal extinction 
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Figure 1. The colour- magnitude diagram of stars in Holmerg II. 
The magenta filled circles at the tip of the main sequence repre- 
sent OB stars used in the analysis. The coloured dots (from blue 
to orange) represent BHeB stars used in this study, colour coded 
to reflect their age (quantified in the colour bar). The red region 
denotes RGB stars used as the background template population 
in the two-point correlation function analysis (see 13. li fe II. It . 



effects. The former primarily influences the colour of the 
locus of BHeB stars, while the latter can introduce broad- 
ening of the BHeB and MS tracks, possibly blending the 
two together in instances of high extinction (e.g., Cannon et 
al. 2003). The galaxies we consider in this paper have been 
selected in part because they have low values of differential 
extinction, which makes for a clear separation between the 
MS and BHeBs. 

The distinct separation between these sequences makes 
it possible to visually select the BHeBs. We select a polygo- 
nal region around the BHeBs, with the blue and red colour 
limits selected to tightly enclose the BHeB sequence (see 
Fig.[T}. The bright end of the BHeBs is determined using the 
10 Myr MS turn off magnitude. For ages < 10 Myr, BHeB 
and MS stars are usually sparse and blend together making 
the separation of sequences not as clear. The faint end of the 
BHeB sequence extends into the bright end of the red clump 
(RC). The precise age and magnitude of where BHeBs in- 
tersect the RC depends upon physical parameters such as 
metallicity and SFH (which can affect the distribution of 
RC stars) as well as the filters used for the observations. 
The faint limit of the BHeB is difficult to identify visually. 
We thus use simulated CMDs to determine the magnitude 
limit at which the BHeBs merge into the RC. Inclusion of 
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some RC stars among our oldest BHeBs stars is still likely. 
We note, however, that RC stars generally have a smooth 
spatial distribution, and thus their inclusion as BHeB stars 
serves to wash out, and not enhance, signals of stellar struc- 
ture. 

Comparison of the observed BHeB sequence with stel- 
lar models allows us to assign ages to individual BHeBs. To 
convert from magnitude to age, we first select the appropri- 
ate metallicity isochrone based on Hn region abundances, 
if available. Otherwise we use the best fit isochronal metal- 
licity as determined by quantitative derivation of the SFH 
(e.g., Dolphin 2002). The stellar models contain the ages of 
the BHeBs sequence tabulated at time intervals of log(age) 
= 0.05. However, our statistical structure analysis is driven 
by having an approximately equal number of stars per time 
bin, and thus we seek a continuous age distribution for in- 
dividual stars, which permits a flexible choice in age bin- 
ning. To create a continuous age distribution, we consider 
the BHeB sequence in log(a<;e)-absolute magnitude space. 
Here, we can functionally approximate the BHeB sequence 
by a 5th order polynomial, which provides for a one-to-one 
analytic relation between age and luminosity. From this re- 
lationship we can convert the obseved magnitudes of BHeBs 
into ages for individual stars. 

Errors on the ages are derived from the photometric 
uncertainties of the individual stars. The bright and faint 
photometric errors for each star are propagated through the 
analytic age-luminosity relationship, and are translated into 
age error estimates. For the galaxies considered in this pa- 
per, the uncertainties in the assigned ages are typically ~ 



10%. This error process does not account for systematic un- 
certainties inherent to the choice of isochrones. Thus, while 
the resultant relative ages are robust, the absolute ages could 
depend on the choice of stellar model. We note that when 
comparing the Geneva and Padova evolutionary models for 
age assignment to BHeBs, Dohm-Palmer et al. (1998) find 
that the two models provide nearly identical results. 

While BHeBs allow us to trace structure from ~ 10- 
500 Myr ago, we can also use young MS stars for additional 
structural information. We select MS stars that are brighter 
than the 10 Myr MS turnoff magnitude and bluer than the 
blue edge of the selected BHeBs. Due to the brightening of 
MS stars during their evolution, some older stars will fall 
into our selected region of colour/magnitude space. From 
Figure 4 in Weisz et al. (2008) we see that stars as old as 
25 Myr may enter our selection, but that they are rare, and 
that the sample is dominated by younger stars. From the 
simulations discussed in § 13. 21 we assign a representative age 
of 8 Myr for stars that pass this selection, and will refer to 
them as OB stars, although we note that our selection does 
not include all true OB stars in each galaxy. 

3.1 Reference Population Selection 

For a comparison reference population, we selected red giant 
branch (RGB) stars from the CMDs. The RGB stars repre- 
sent an older stellar component that has had sufficient time 
to diffuse and erase signs of local structure. 

The reference population was selected by drawing a box 
in the RGB at least 0.5 mags fainter than the tip of the red 
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giant branch (TRGB). The colour limits were restricted to 
be redder than the RC, to avoid possible selection of RHeBs, 
slightly bluer than the red edge of the RGB. An example of 
the RGB stars selected in Holmerg II is show in Fig. [T] 

3.2 An Example of the Star Selection: Holmerg II 

In this section, we describe the star selection method 
for the CMD of the M81 Group dwarf irregular galaxy 
Holmerg II. The CMD presented in Fig. [T] shows the MS, 
BHeB, and RGB stars used for tests of structure evolution 
in Holmerg II. 

The selection of BHeBs requires setting both colour and 
magnitude limits around the BHeB sequence. Holmerg II, 
like many low metallicity galaxies, shows a clear separa- 
tion between the BHeBs and MS stars over several mag- 
nitudes (rriF555w ~ 21-25). For slightly fainter magnitudes, 
the MS broadens due to photometric effects, and the BHeB 
sequence curves sightly redward. However, it is possible to 
visually identify the separate sequences and assign a blue 
colour limit. The red colour limit is also determined by eye 
and is based upon the decreasing density of BHeBs as a 
function of increasing colour. The luminous stars between 
the BHeB and red helium burning (RHeB) sequences are 
primarily HeBs transitioning between the two phases. Lu- 
minosity evolution during this transition is roughly constant 
(Marigo et al. 2008), thus including them does not strongly 
affect the results. Main sequence turn-off (MSTO) stars are 
known to also cross the BHeB sequence, however, the rela- 
tive number of stars is only a few percent (Dohm-Palmer et 
al. 1997), and does not significantly affect the BHeB selec- 
tion. 

The BHeBs merge into the RC at rriFhhhw ~ 26.5. In- 
clusion of some RC stars in the BHeB selection is proba- 
ble, due to confusion at the faint limit. However, RC stars 
are several Gyr in age, and do not exhibit clustered spatial 
structure. Thus, if counted as BHeBs, they tend to smear 
out structure, rather than enhancing it. To quantify the faint 
BHeB limit, we constructed a simulated CMD for a constant 
SFH, which was convolved with the Holmerg II artificial star 
tests to account for observational and completeness effects. 
We then applied the same colour selection selection criteria 
to the simulated CMD, for which we know all stellar ages 
and masses, and determined the spread in BHeB age as a 
function of m,F555W- Typically, the age spread at a given 
magnitude for BHeBs is ~ 10%. However, at ttifsssw ~ 27 
(~ 550 Myr), the spread increases to 20% and continues in- 
creasing thereafter. Although we have selected BHeBs stars 
in Fig. [l]older than this limit, we note the uncertainties at 
the oldest ages are as high as 40%. 

BHeBs are only reliable age indicators for ages > 10 
Myr. To gain additional information on structure in the 
most recent 10 Myr, we selected luminous MS stars using 
the isochrones of Marigo et al. (2008). The faint limit is 
based upon the 10 Myr MSTO isochrones, Mfsssw ~ —4. 
However, MS luminosity evolution causes slightly older stars 
to become brighter than this limit. Although this is a rapid 
phase in a star's evolution, we chose a more conservative 
values of Mfsssw = —5 (rriFsssw = 22.65 - at the average 
distance to our targets) to minimise sources of contamina- 
tion. The red colour limit of the OB stars was set by the 
blue edge of the selected BHeBs. We also required MS stars 



to have a colour of rriFsssw — rriF8i4w < 0.5. There are ~ 
10 objects blueward of this limit, which could be binaries 
systems or unresolved blends, however they are not likely 
single OB type stars. 

We have also carried out consistency tests to check 
the number of observed vs. expected OB stars. We simu- 
lated a CMD (30 realisations) with the star-formation his- 
tory of Holmerg II based on the analysis of W08 for the 
past 10 Myr. A SFR of 0.17 M yr _1 was adopted for the 
past 4 Myr based on the Ha flux and the SFR from 4 to 
10 Myr was derived from the analysis of the BHeB stars 
(see W08). We adopt a stellar IMF with index, -1.30 (i.e. 
close to the Salpeter value), from 0.15 to 100 Mq that cov- 
ers the full range of masses for 10 Myr. This yields 8.8 * 10 5 
stars, 160± 12 of which pass our colour/magnitude selection 
criteria of being a high-mass MS star. The lowest mass star 
that passes our criteria has a present day mass of ~ 20 Mq . 
We observe 126 stars in our selection box indicating that 
not all high-mass stars in the galaxy made it through our 
selection process. However, we only accepted 'point-like' ob- 
jects in our analysis, so all stars in clusters or associations 
would have been removed from the sample. Additionally, 
some stars may have been removed from the sample due 
to the complicated background found in Hn regions. BHeB 
stars, being older and more distributed on average, are not 
expected to suffer from these selection effects. We conclude 
that we are picking up the majority of the expected MS mas- 
sive stars. The mean and median age of the stars meeting 
our criteria in selecting massive MS stars is ~ 7 — 8 Myr. 
We therefore adopt this age as a representative value for our 
MS population. 

To provide a reference population, we sub-sampled the 
RGB using a simple polygon. The bright and faint limits 
were selected to run parallel to the TRGB, while the red and 
blue limits were chosen to to minimise RHeB, Red Clump 
(RC), and Asymptotic Giant bBranch (AGB) star contam- 
ination. 

The CMD in Fig. Q] show the OB stars (purple points; 
126 stars), BHeBs (age coloured sequence; 23,653 stars), and 
RGB reference population (red points; 28,919 stars). The 
corresponding spatial distribution of these populations are 
shown in Figure(3] The position and derived age of all BHeBs 
and young stars is given in Table [T] The error in age for each 
star is derived from its photometric error and does include 
uncertainties in the models used. 



4 STRUCTURAL ANALYSIS 

Once the age of each star has been estimated, we can begin 
to study their spatial distribution as a function of age. For 
this, we break the sample up into discrete age bins, and 
apply statistical tests to each distribution. We apply two 
tests to the data, each will be discussed in some detail below; 
1) the two-point correlation function (TPCF) and 2) the Q- 
parameter (Cartwright & Whitworth 2004). 

In Fig. [3] we show an example from the survey, 
Holmerg II. In each panel we show the spatial distribution 
of 500 stars (except for the first panel where only 126 stars 
are used) stochastically sampled with the mean age given for 
four age bins selected from the BHeB sample, as well as for 
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Table 1. The position, age and estimated age error for all 
BHeB and young stars in our sample. Stars that have a 'MS' 
after their ID number were selected as young main sequence 
stars (see §|3}. The full version is available in the electronic 
version of the paper. 



Galaxy 




star ID 


X 

(pixels) 


y 

(pixels) 


age 
(Myr) 


age error 
(Myr) 


Holmbcrg 


II 


1 


2052.95 


2274.90 


6.7 


0.2 


Holmbcrg II 


2 


2164.82 


1628.10 


7.2 


0.2 


Holmbcrg 


II 


3 


1459.65 


2061.96 


6.7 


0.2 


Holmbcrg 


II 


4 


2368.83 


1843.51 


7.7 


0.2 


Holmbcrg 


II 


5 


2858.23 


2619.56 


8.4 


0.2 


Holmbcrg 


II 


6 


2473.65 


2227.54 


8.1 


0.2 


Holmbcrg 


II 


7 


2750.72 


2117.37 


9.0 


0.2 


Holmbcrg 


II 


8 


2639.68 


624.23 


9.1 


0.2 


Holmbcrg 


II 


9 


1625.98 


1657.29 


9.4 


0.2 


Holmbcrg 


II 


10 


2752.44 


1774.87 


9.6 


0.2 



the OB stars and RGB star (reference) sample. Examples 
for each galaxy are given in Appendix El 

4.1 Two point correlation function (TPCF) 

The TPCF is a relatively straight forward way to study 
spatial structures. This method determines the distance be- 
tween all possible pairs of stars (d s ), shown as a histogram, 
which is then normalised to that of a reference distribution, 
i.e., Nii nks (d s )/N re ference(ds). Gomez et al. (1993) used the 
TPCF to study the distribution of young stars in a star- 
forming region and used a random distribution for the refer- 
ence distribution. This is not applicable in the current case 
as the galaxies have underlying structure, being centrally 
concentrated. For the reference distribution we use an equal 
number of stars drawn from the red-giant branch, where 
we can be confident that each star has an age greater than 
~ 1 Gyr. 

For the TPCF, an equal number of stars per age bin 
(and also for the reference distribution) is required. For the 
galaxies presented here, we use between 500 and 1250 stars 
per age bin, depending on the total number of BHeBs stars 
present. The exact number was chosen for each galaxy as a 
compromise in age resolution and having enough stars for 
the statistical analysis. For all galaxies, multiple numbers of 
stars per age bin were used in order to check the robustness 
of the results. For each age bin, we generated five reference 
distributions, randomly sampling the same number of stars 
from the background population. The TPCF was then de- 
rived with each of the five reference distributions and the 
slope and y-offset of the resultant histogram was measured. 
The adopted values were the mean of the five distributions 
with the standard deviation being used to estimate the un- 
certainties. 

Ideally, the slope and y-offset of the TPCF should con- 
verge to a value of zero if the population under study is 
evolving into the reference distribution. Deviations from zero 
imply that the reference population is not an ideal fit, where 
the offset is due to contamination or due to a real offset. 
Differing distributions of young and old stars have been re- 
ported in the literature before (e.g., UGC 4438 - Dolphin 
et al. 2001). However, since we will be studying relative 
changes, the adopted reference population is not critical. 



Since young stars are known to be more highly clus- 
tered relative to older stars, we expect the TPCF to have 
the largest offset from the reference population at young 
ages and evolve towards the reference distribution as older 
populations are considered. Ideally, we should see both the 
slope and the y-offset of the TPCF stop evolving at some 
characteristic time, t CV o, where all of the original structure 
has been erased on galactic scales. The measured TPCF for 
each galaxy in the sample is shown in the middle and bottom 
panels of the evolutionary figures (Figs. Q] [5] [6] [3 [H & 0. 
Finally, in order to test for binning artefacts, we have also 
adopted smaller numbers of stars included in each age bin. 
These bins, which have increased noise but higher time res- 
olution, are shown as light filled squares in the evolutionary 
panels. 

4.2 Q-parameter 

Cartwright & Whitworth (2004; hereafter CW04) have de- 
veloped a technique to measure the degree of substruc- 
ture/degree of central concentration of a discrete two- 
dimensional dataset. We refer to the reader to Cartwright & 
Whitworth (2004, 2009), Schmeja & Klessen (2006), and B09 
for a full description and tests of the technique. The CW04 
method, known as the Q parameter, is based on Minimum 
Spanning Trees (MST), which for each dataset is formed by 
connecting all points (spatial positions in this case) in or- 
der to form a unified network, such that the total length 
(i.e., sum) of all of the connections, known as 'edges' or 
'branches', is minimised, and no closed loops are formed. 
The mean branch length of the MST is calculated and then 
normalised by the factor (-/Vtotaij4)°' 5 /(A?totai — 1) to give m, 
where A is the area of the dataset and iV tota i is the total 
number of stars in the dataset. The mean distance between 
all points (derived from the two-point correlation function) 
is then computed, which is then normalised by the radius of 
the dataset, giving s. In our case, the area is derived from 
the convex hull of the dataset, and the radius is assumed to 
be the radius corresponding to a circle with that area (see 
Schmeja & Klessen 2006). Q is then simply fh/s. 

CW04 have carried out a series of numerical tests and 
have shown that Q is highly sensitive to the amount of struc- 
ture present (i.e., the degree of substructure) in a dataset. 
Additionally, Q can differentiate between a highly structured 
and a centrally concentrated distribution. Their tests used 
fractal distributions to simulate clustering, where low frac- 
tal dimensions have large amounts of substructure, and high 
fractal dimensions are more uniformly distributed. CW04 
found that a 3D spherical distribution (projected onto a 2D 
surface) that is uniformly populated will have Q = 0.79, 
while highly structured distributions will have Q < 0.79 and 
centrally concentrated distributions will have Q > 0.79. 

In general, the lower the Q-value, the more substruc- 
ture is present. As the dataset becomes more uniform or 
centrally concentrated, Q increases. One complication is the 
effect of elongation on the value of Q. B09 and Cartwright 
& Whitworth (2009) have shown that if a dataset is elon- 
gated, the Q value will be lower than for the same spherical 
distribution. This is important for the current work, given 
that many of the dlrrs in the current study are highly elon- 
gated and are seen from non-face-on orientations. However, 
in the current work we shall not use the absolute value of 
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the measured Q-value, but instead study the evolution of Q 
as a function of age, which should be largely independent of 
elongation and other geometrical effects. 

Additional tests on the effects of differential extinction 
and contamination on the measured Q-value have been dis- 
cussed in detail in Appendix B and C of B09, respectively. 
These are not expected to significantly affect the results pre- 
sented in the current study. We will focus our attention on 
relative changes in Q, hence (except for the youngest stellar 
populations) all populations should be similarly affected by 
extinction throughout the galaxy. 

As with the TPCF, due to the highly structured distri- 
bution of young stars in galaxies, we expect Q to start with 
a low value and evolve towards the reference distribution 
and stop evolving once t eV o is reached. i ovo was estimated 
for each galaxy using the evolutionary figures (Figs. [4] [5j [6j 
[7] [8] & [9} where the distributions stop changing. This was 
done by eye, weighting the three measurements. 

Tests of the stability of Q using the observed distribu- 
tion as well as models showed that fewer sources are needed 
than for the TPCF method. An additional benefit of the Q- 
parameter method is that it is independent of the number 
of sources used. We take advantage of this feature by using 
smaller bins at young ages (where we expect the fastest and 
most pronounced structural evolution) and larger bins at 
older ages. Typically, between 500 and 750 stars were used 
in the first five age bins, and 1500 stars in the oldest bins. 
These points are shown as filled circles in the top panels of 
the evolutionary figures (Figs. H [S] H H & 0. Since a 
constant number of sources is not necessary, we have also 
measured the Q-value for the OB stars in our sample. This 
is shown as a filled triangle in the evolutionary figures. 

In order to test the reliability of the results, and any 
dependence on binning, we have also carried out the same 
analysis using smaller age bins. While this increases the 
noise in the results, it can be used to confirm general trends. 
These points are shown as light-shaded squares in the evolu- 
tionary figures. Finally, we have estimated the uncertainty 
associated with the determination of the Q-value by pick- 
ing N stars at random from the reference population of 
Holmerg II. We find that the uncertainty in Q is 0.027, 
0.014, 0.011, 0.010 for TV = 200, 500, 1000, & 1500, respec- 
tively (based on 100 random realisations of the background) , 
which are shown in the figures as error bars. 



4.3 Caveats 

While the two methods employed in the current work are 
powerful tracers of structural evolution, some caveats are 
necessary to discuss. We have already discussed the need for 
using a representative reference distribution for the TPCF 
method (and to a smaller degree, the Q- method). Below are 
a few additional caveats. 

(i) The existence of long lived structures (i.e., stel- 
lar clusters): Not all structures are expected to be de- 
stroyed within the first few hundred Myr of their lives. Stel- 
lar clusters are gravitationally bound collections of stars that 
may live to very old ages. However, only a very small frac- 
tion of stars in a galaxy are in clusters (after a few Myr - e.g., 
Miller & Scalo 1978; Adams & Myers 2000; Goddard, Bas- 



tian, & Kennicutt 2010). Hence, stellar clusters are unlikely 
to significantly affect the results presented here. 

(ii) Inclination effects: The galaxies in our sample are 
observed at various inclination angles. While Q has been 
developed in order to give a quantitative measure of the 
amount of structure present, the inclination angle and elon- 
gation of the system complicate the interpretation. Hence, 
instead of absolute values, we will focus on relative changes 
within single galaxies. In this way we can track the evolution 
of structures, but cannot quantify the fractal dimension of a 
population or the level of central concentration for the older 
populations. 

(iii) Limited field of view: In addition to inclination 
effects, some galaxies used in our survey have stellar pop- 
ulations that continue beyond the ACS coverage. This also 
limits the use of the Q-value as an absolute measure. How- 
ever, the evolving spatial structure should be evident within 
the observed field-of-view of each galaxy, again meaning that 
we can track the relative change of structure. 



5 RESULTS 
5.1 Holmberg II 

Holmberg II is an actively star-forming dwarf galaxy in the 
M81 group and has an inclination angle of i ~ 44 degrees 
(Puche et al. 1992). Its current SFR (SFR~ 0.17 Meyr -1 ) 
is roughly a factor of two higher than the average over the 
previous 500 Myr (W08), suggesting a recent increase in its 
activity. The HI distribution within Holmerg II is made up 
of numerous holes, although the distribution of red (presum- 
ably older) stars within the galaxy is quite smooth (W08). 
The ACS observations used here do not cover the full ex- 
tent of the galaxy, as both red and blue stars extend off 
the northwest of the field-of-view. Due to the recent active 
star-formation within this galaxy, numerous high-mass (OB) 
stars are seen on the main-sequence (see Fig. [1} . 

Figure [4] shows the evolution of the Q-value, TPCF 
slope, and TPCF y-offset for the BHeB stars (filled circles) 
and OB stars (filled triangle - only in the Q-value plot). Note 
the strong evolution in all indicators until an age of 225 Myr. 
The distributions then level off, and remain approximately 
constant until age of ~ 600 Myr (our last age bin). 

One important point to note is that the Q-value does 
not merge into the reference population nor does the TPCF 
converge to zero. The reason for this is that the reference 
population (RGB stars) has a significantly different distri- 
bution than the BHeB stars at all ages. This holds for all 
galaxies in our sample, meaning that we will limit our anal- 
ysis to the relative change in each of the structural analysis 
methods. In the case of Holmerg II the RGB stars appear to 
be more centrally concentrated than the BHeB stars, which 
is opposite to what we see in all of the other galaxies in our 
sample. 

In Fig. [3]we show the spatial distribution of BHeB stars 
in Holmerg II for six age bins. The strong evolution in the 
clustering properties of the stars can be discerned by eye, 
which is quantified in Fig. [4] 
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Figure 3. The positions of the BHeB stars in Holmerg II for five age bins and a random selection from the reference RGB stars. Each 
panel has 500 stars (except for the first, where we only have 126 OB stars) that were stochastically chosen from a given bin. The mean 
age of the bin and the measured Q-value is given in each panel. 



5.2 DDO 165 

DDO 165, a dwarf galaxy in the M81 group, was covered 
by a single ACS pointing. The CMD of DDO 165 shown in 
Fig.[2]does not show a clear separation between BHeB stars 
and the main- sequence, showing just a single blue plume. 
However, given the post-burst nature of this galaxy (Lee et 
al. 2007, W08), the MS is poorly populated at the brightest 
magnitudes compared to the BHeBs. Brighter blue stars are 
more likely BHeBs, making identification trivial. However, 
at fainter magnitudes, separation of the MS and BHeBs be- 
comes more difficult. The evolution of structure as traced by 
the BHeB stars can be seen in Fig. [5] As with Holmerg II, 
we see a clear evolution in structure at young ages with all 
indicators. We adopt a value of 350 ± 75 Myr for t cvo for 
DDO 165, although we note that the observations are also 
consistent with a continuously evolving population until the 
oldest age bin (~ 500 Myr). 

5.3 IC 2574 

IC 2574 is a large, highly elongated, dwarf Irregular galaxy 
in the M81 group (i ~ 53 degrees - Oh et al. 2008). The 
galaxy has a ring of bright Hn regions around an Hi hole in 
the northeast of the galaxy, the brightest of which contains 
~ 1/3 of the ionising luminosity as 30 Doradus (Walter et 
al. 1998). Three ACS fields were used to map its structure. 
As seen in W08, the stellar distribution within IC 2574 is 



highly non-symmetric within the field-of-view, especially at 
young ages (blue stars). These are concentrated largely in 
the northern-most ACS field, while the red stars are more 
symmetrically distributed (being most highly concentrated 
in the centre field-of-view). IC 2574 appears to have under- 
gone a period of intense star-formation, beginning approxi- 
mately 1 Gyr ago and lasting for ~ 500 Myr (W08). 

The evolution of the spatial structure in IC 2574 is 
shown in Fig. [U] In the upper panel, Q evolves gradually 
from the youngest age bin (the OB stars) until ~ 150 Myr, 
at which point it decreases until an age of ~ 200 Myr at 
which point it begins to rise again until the oldest age bin. 
We see the same type of evolution in the TPCF method, 
although it is less apparent at young ages due to the higher 
number of sources needed in order to carry out the analysis 
(i.e., causing the age resolution to decrease). We have used 
different age binning choices for the TPCF and find that the 
results corroborate the findings from the Q evolution. 

We conservatively adopt t cvo = 150 Myr, although the 
structural properties of IC 2574 clearly continue to change 
until the oldest ages that we can study (~ 500 Myr). Hence 
our value should be taken as a lower limit. 

5.4 NGC 2366 

NGC 2366 is the final dwarf galaxy in our sample from the 
M81 group. Two ACS fields were required to observe the 
main optical extent of the galaxy. The older stars show an el- 
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Figure 4. The evolution of spatial structures in Holmerg II. Top 
panel: The measured Q-value as a function of age. Filled cir- 
cles represent the mean age and Q-value of BHeB stars while 
the filled triangle is that for OB stars selected from the CMD. 
Lighter square symbols denote the results using more refined age 
bins, i.e. including fewer stars per bin. The vertical dashed line 
is the adopted age (with uncertainty indicated as the filled area) 
where structural evolution stops. The solid and dashed horizon- 
tal lines represent the value and error of the RGB background, 
respectively. Middle panel: The same as the top panel but now 
we show that evolution of the slope of the two-point correlation 
function (TPCF). Bottom panel: The same as the middle panel, 
but now the y-offset of the TPCF. Note the strong evolution in 
Q and the TPCF until ~ 225 Myr, followed by a flat distribution 
until an age of ~ 600 Myr. 



liptical, centrally concentrated structure that nearly fills the 
field-of-view, while the younger (blue) stars show a highly 
structured distribution (W08). The galaxy contains a large 
Hn region, NGC 2363 (Kennicutt et al. 1980) which con- 
tains a large population of young massive stars. NGC 2366 
appears to be undergoing a recent increase in its SFR, which 
began ~ 50 Myr ago (W08). The inclination angle of the 
galaxy is i ~ 64 degrees (Oh et al. 2008). 

In Fig. we show the evolution of spatial structure 
within NGC 2366. The large amounts of structure seen at 
young ages rapidly disappears. From the evolution of Q, 
we see that the distribution becomes flat (i.e. lacks further 
change) at ~ 100 Myr. The TPCF slope and y-offset suggest 
a more complicated picture, changing until ~ 200 Myr and 
then reversing the trend. We do not have an explanation for 
this behaviour, and we adopt t cvo from the Q evolution as 
a conservative lower limit. 



Figure 5. The same as Fig. H] but now for DDO 165. 
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Figure 7. The same as Fig.|3]but now for NGC 2366. 



Figure 8. The same as Fig. 2]but now for NGC 784. 



5.5 NGC 784 

NGC 784 is a barred dwarf Irregular galaxy seen nearly edge 
on. The galaxy hosts a number of Hll regions showing that 
star-formation is still actively ongoing (Kandalyan, Khas- 
sawneh, & Kalloghlian 2003). Due to its highly inclined po- 
sition (i ~ 79 degrees, Karachentsev et al. 2004), it is more 
difficult to see the evolution of structures by eye, although 
the TPCF and Q evolution suggest that it is happening. 

The TPCF and Q evolution of BHeB stars, shown in 
Fig. [S] shows a smooth evolution from a highly structured 
distribution to a more uniform/centrally concentrated type 
grouping. There is some evidence that the distribution stops 
evolving at an age of 225 Myr, although it may continue past 
the last age bin available to us. The young OB stars are 
consistent with the extrapolation of the Q-parameter from 
older ages. 



5.6 NGC 4068 

NGC 4068 is a blue compact dwarf (BCD) galaxy located 
at a distance of 4.2 Mpc, in the Canis Venatici I (CVn I) 
cloud (Karachentsev et al. 2006). The galaxy is currently 
undergoing an increase in its star-formation activity, which 
began ~ 400 Myr ago and is continuing until the present 
(McQuinn et al. 2009). NGC 4068 has an inclination angle 
of ~ 60 degrees (Karachentsev et al. 2004). 

Figure [9] shows the evolution of stellar structures in 
NGC 4068. While the TPCF measurements are quite noisy, 
the Q evolution appears to be quite smooth. We assign a 
value of t cvo > 325 ± 50 Myr, although there is not strong 



evidence that the evolution of structure stops at this age. 
While there are only a small amount of OB stars appar- 
ent in this galaxy (60), their position in the Q-diagram (top 
panel Fig. O is consistent with the extrapolation of the Q 
evolution seen in the BHeBs. 



6 DISCUSSION AND CONCLUSIONS 

Each galaxy in the present work shows clear evidence of 
evolving stellar structures. Stars appear to be born in a 
highly structured, possibly hierarchical or fractal distribu- 
tion, which evolves smoothly into a more uniform (centrally 
concentrated but lacking significant substructure) distribu- 
tion. Using two statistical methods, the Q parameter and 
the Two Point Correlation Function (TPCF), we have been 
able to place lower limits on the timescale over which this 
transition occurs. We find that this transition takes between 
90 Myr (NGC 2366) to > 300 Myr (DDO 165). These are 
similar timescales to those that we found in previous stud- 
ies on the LMC (175 Myr - B09) and SMC (75 Myr - G08). 
Hence, stellar structures remain coherent within these galax- 
ies, on average, for over 200 Myr. 

The estimated evolutionary timescales, t ovo , are lower 
limits as we see the TPCF and Q-parameter continue to 
vary to some degree in all galaxies (with the exception of 
Holmerg II) until the oldest age bin. Some galaxies, such as 
DDO 165, NGC 784, and NGC 4068 appear to have contin- 
uously evolving structures until the oldest ages under study. 
Others, such as NGC 2366 and IC 2574 appear to show 
more of a composite evolution, with a rapid initial evolution 
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Figure 9. The same as Fig.|4]but now for NGC 4068. 



(100 — 150 Myr), followed by either a period of no or oppo- 
site evolution, and finally a further decrease in substructure 
at large ages. Holmerg II shows the cleanest evolution with 
a rapid decrease in structure lasting ~ 225 Myr with no fur- 
ther evolution found until the oldest age bin in the sample 
(~ 600 Myr). 

There does not appear to be a strong dependence on 
galaxy luminosity or galaxy size. Unfortunately, stellar ve- 
locity dispersion measurements are not available for these 
galaxies. For the LMC and SMC, where such measurements 
are available, B09 and G08 found that t ovo was comparable 
to the crossing time of the galaxy. Such a prediction can be 
verified with the current study once the velocity dispersion 
of the young stellar components are measured. Using the 
measured physical size of the galaxies and adopting the de- 
rived evolutionary timescales, we can estimate upper limits 
to the stellar velocity dispersions of these galaxies, assum- 
ing that structures erase on a crossing time. As an example, 
our field of view of Holmerg II is approximately 360" by 
200". Approximating this area as a circle leads to a radius 
of ~ 150" (or 2.46 kpc assuming a distance of 3.39 Mpc - 
Karachentsev et al. 2002). Substituting t cvo for the crossing 
time (icross = R/a) leads to the upper limit of the velocity 
dispersion of ~ 11 kms -1 . As a substitute for the stellar ve- 
locity dispersion, we note that the mass-weighted median of 
the HI velocity dispersion in Holmerg II is 9.4 km s _1 (Tam- 
burro et al. 2009). This may be taken as an upper bound 
for the velocity dispersion of the newly formed stars, as the 
measured velocity dispersions show both cold (narrow) and 
broad (warm) components (see Young and Lo 1997 and dis- 
cussion in Tamburro et al. 2009), and the star formation is 



expected to be associated with the narrow component. Thus, 
for Holmerg II there appears to be no discrepancy with the 
velocity dispersion measurements and the existence of sub- 
structures corresponding to the crossing time. For smaller 
galaxies, the crossing time argument becomes more restric- 
tive, but, with these observations, we do not see a trend in 
structure lifetime with galaxy size. 

In nearly all cases, the RGB distribution (stars with 
ages >1 Gyr) is significantly different than our oldest age bin 
accessible with the BHeB method, with the RGB stars be- 
ing more extended. This is consistent with previous findings, 
such as UGC 4438 (Dolphin et al. 2001) and NGC 6822 (de 
Blok & Walter 2003) and similar to what we found for the 
LMC (B09) and SMC (G08). Population studies and numeri- 
cal modelling of galaxy formation and evolution (see Stinson 
et al. 2009) have shown that the outer parts of dlrrs are older 
on average than the inner parts. At older ages, part of the 
evolution that we have detected in our sample of dlrrs may 
be due to this effect, and may not be directly attributable to 
dynamical structural evolution. The one counter-example in 
our sample is Holmerg II, where the RGB stars appear to be 
more centrally concentrated than the BHeB stars. Without 
considerably increasing the field-of-view of our observations, 
it is difficult to speculate why the distribution of RGB stars 
(relative to the younger BHeB stars) is different than for the 
other galaxies in our sample. 

The main limitation of the current analysis is the lack 
of information on how stellar structures are dissolving on 
different spatial scales, as we have limited our study to the 
global distribution in each galaxy. Studies that are able to 
probe the evolution of structure in both dimensions (in time 
and as a function of spatial scale) will be able to definitively 
test what mechanism is driving the observed evolution. 
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APPENDIX A: SPATIAL EVOLUTION OF 
STELLAR STRUCTURES FOR EACH GALAXY 
CONSIDERED 

Here we present examples of the spatial evolution of stel- 
lar structures for each galaxy in our sample. Each figure is 
similar to Fig. [3] In all cases, younger populations are seen 
to have much more substructure which dissolves in tens or 
hundreds of Myr. 

This paper has been typeset from a TjtX/ IMeX file prepared 
by the author. 
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Figure Al. The positions of the BHcB stars in DDO 165 (similar to Fig. [3jl for five age bins and a random selection from the reference 
RGB stars. Each panel has 500 stars that were stochastically chosen from a given bin. The mean age of the bin and the measured Q- value 
is given in each panel. 
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Figure A2. Similar to Fig. [3] but now for IC 2547. 
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Figure A3. Similar to Fig. [3] but now for NGC 2366. Only 300 stars per age bin are shown in order to highlight the structural evolution 
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Figure A4. Similar to Fig. [3] but now for NGC 784. Only 400 stars per age bin are shown in order to highlight the structural evolution 
at early times. Note that the edge on structure of the galaxy is likely to washout the the signal of structural evolution. 



© 2009 RAS, MNRAS 000.ITHT21 



The evolution of stellar structures 15 



4000 
3000 
.* 2000 

Cl 

1000 



■ 49 Myr £ 

Q = 0,700 NGC 4068 


' i " 1 ■ 1 " 1 1 

' 90 Myr. .. 

• • • Mif&f 

= 0.772 


■ i ' 1 ■ n ■ i 1 

126 Myr t f « 

Q = 0.821 


■ 1 1 1 — tH 1 1 1 

' 184 Myr 

* * ** * 
■ ■ . -* * 

Q = 0,832 

■ ■ 


— ■ — i — i — i — i — i — ■ — i — i — 
' 234 Myr , . . 

Q = 855 


'RGB stars'-* *% " 

ft* • V ' * 
* . .. * !• '.v ■*■ . 

■ ; *.-..» * : 

■ * 1 "**t ■ *. 

q =ue5&.' *■ 



2000 4000 

pixels 

Figure A5. Similar to Fig- [3j but now for NGC 4068. Only 400 stars per age bin are shown in order to highlight the structural evolution 
at early times. 
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